Bias asymmetry in the conductance profile of magnetic ions on surfaces probed by 

scanning tunneling microscopy 
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The conductance profiles of magnetic transition metal atoms, such as Fe, Co and Mn, deposited 
on surfaces and probed by a scanning tunneling microscope (STM) , provide detailed information on 
the magnetic excitations of such nano-magnets. In general the profiles are symmetric with respect to 
the applied bias. However a set of recent experiments has shown evidence for inherent asymmetries 
when either a normal or a spin-polarized STM tip is used. In order to explain such asymmetries here 
we expand our previously developed perturbative approach to electron-spin scattering to the spin- 
polarized case and to the inclusion of out of equilibrium spin populations. In the case of a magnetic 
STM tip we demonstrate that the asymmetries are driven by the non-equilibrium occupation of the 
various atomic spin-levels, an effect that reminds closely that electron spin-transfer. In contrast 
when the tip is not spin-polarized such non-equilibrium population cannot be build up. In this 
circumstance we propose that the asymmetry simply originates from the transition metal ion density 
of state, which is included here as a non-vanishing real component to the spin-scattering self-energy. 



I. INTRODUCTION 

The possibility of altering and controlling the spin- 
state of a single magnetic ion or of a small magnetic 
cluster with an external probe represents a unique oppor- 
tunity towards the understanding and the exploitation of 
the magnetic interaction at the most microscopic level. 
Possible areas of application for such ability may include 
spin-based quantum logic, where one necessitates to pre- 
pare, manipulate and read spin-qubits. It is then crucial 
to develop tools capable of addressing the single spin- 
limit. Low-temperature scanning tunneling microscopy 
provides one of such tools. In general the method ex- 
ploits a scanning tunneling microscope (STM) operated 
in spectroscopical mode, by which the inelastic electron 
tunneling spectroscopy (lETS) at the spin-excitations of 
a given system is measuredM^. This scheme is known as 
spin-flip lETS (SF-IETS). The same STM can also be 
used to position and manipulate the magnetic atoms on 
a non-magnetic substrate^, so that STM appears both as 
a fabrication and subtle characterization tool. 

Transition metal magn eti c a toms on insulating sur- 
faces, in particular MiPl, CcP^and F^^, have been the 
focus of intensive research in the last few years. These 
have all partially filled d-shells, which are highly localized 
and responsible for the magnetic moment, and extended 
s-like electrons, which are responsible for the electron 
conduction. In general s and d electrons interact via ex- 
change coupling so that the magnetic structure is coupled 
to the conducting electrons. The magnetic atoms are usu- 
ally deposited on carefully prepared CuN-decorated Cu 
surfaces, where the typical electronic coupling is weak 
enough that the magnetism is preserved, but it is suffi- 
ciently strong to break the atomic central symmetry so 
that magnetic anisotropy develops. STM experiments are 
then conducted and the fingerprint of a magnetic excita- 
tion is a step in the differential conductance, G — dl/dV, 
as a function of bias, V {I is the STM current). These 
appear at the critical voltage necessary to open a new in- 



elastic transport channel, i.e. at voltages corresponding 
to the given magnetic excitation energy. 

Several methods aimed at modeling SF-IETS have 
been recently developed. Early theoretical work has fo- 
cused on second order perturbation theory to describe the 
experimental conductance spectra of equilibrium spins 
by either using a master equation approaclpSHlSl or a 
non-equilibrium Green's function one^^. More recently 
this scheme has been extended to third order, which al- 
lows us to describe additional features in the G{V) line- 
shape that cannot be accounted for at the second order 
levep2"22l. These works have been very successful in de- 
scribing conductance profiles, which appear symmetric 
with the external bias polarity, i.e. that they cannot 
distinguish whether the current flows from the sample 
to the tip or in the opposite direction. However, recent 
experiments have shown that regardless on whether a 
non-spin-polarizecP or a spin-polarizecP'^l tip is used the 
lETS profiles exhibit an intrinsic asymmetry with respect 
to the applied voltage, i.e. G{V) ^ G{—V). 

In the case of a spin-polarized STM tip, where the 
tip density of states is spin split between majority (spin 
up) and minority (spin down) carr iers, t he asymmetry 
has been theoretically well explainecpI^EilMl^ j^ j^^^ been 
shown that spin selection rules enforce a suppression of 
the inelastic scattering, which depends on the direction 
of the electrons flow. This results in a asymmetric con- 
ductance profile, where the magnitude of the asymme- 
try depends directly on the spin-polarization of the tip. 
It is also well understood that by driving spins out of 
equilibrium (e.g. by decreasing the tip-sample distance) 
the conductance line-shape changes^EElEMSI, jn this 
case we must assume that the tunneling electrons in- 
fluence the spin-state of the atom as the time between 
inelastic events is small compared to the spin relaxation 
time. A tunneling electron can then encounter the lo- 
cal spin in an excited state far from the ground state. 
The non-equilibrium population of the various accessible 
spin-states then becomes bias-dependent and, for spin- 



polarized tips, this enhances the asymmetry of the G{V) 
hne-shape. 

Also in the case of a non-spin-polarized tip a bias asym- 
metry has been revealed experiment all ji^. In particular 
this appears to be quite prominent for both single Mn 
atoms and Mn mono-atomic chains. This feature has 
been previously ascribed to a shift in the magnetic atom 
on-site energy, i.e. to an effect arising from the details 
of the density of states of the atom producing scattering. 
Such a density of state effect produces a non-trivial slope 
in the conductance as a function of biaa^Sl. The on-site 
energy shift however does not account for the asymme- 
try seen in the inelastic step heights, which also depends 
on bias. Here we provide an alternative theoretical de- 
scription, which allows us to better fit the experimentally 
found conductance line-shape. 



In our previous worksISIiZl we have com bined the 
non-equilibrium Green's function (NEGFpHSSI formalism 
with a perturbative expansion of the electron-spin inter- 
action in order to describe SF-IETS spectra in a manner, 
which is fully amenable to a n imp lementation within den- 
sity functional theory (DFT^^^ESl. The scheme essentially 
consists in constructing an electron-spin interaction self- 
energy, which describes the inelastic tunneling events. 
The interacting self-energy was previously expanded first 
up to second ordei*!^ and then to the third ordepl^, with 
this latter describing the logarithmic decays of G{V) at 
each conductance step. Although both non-equilibrium 
effects and spin-polarized lETS have been well described 
up to the second order by the master equation approach, 
in this work we extend our formalism to include the extra 
line-shape features that can be ascribed to the third or- 
der self-energy. We also propose that the real part of the 
interacting self-energy, which has be en we ll studied in the 
case of electron-phonon interactiona^iESl^ jg g^ necessary 
addition to the NEGF formalism in order to account for 
asymmetric features in non-spin polarized systems. 

The layout of the paper is as follows. In the next 
section we extend our NEGF formalism to account for 
spin polarized leads. In the same section we also derive 
a second-order electron-spin self-energy, which includes 
a second order expansion of the spin propagator. This 
provides the means to study non-equilibrium effects that 
result from high current densities. Furthermore, in the 
case of non-spin polarized tips we derive an expression 
for the real part of the scattering self-energy up to sec- 
ond order. All of the above is combined with the third 
order contribution to the electron propagator. Then we 
move to the results. First we study the non-equilibrium 
effects arising from an increase in current density for the 
case of a Mn dimer. Then we recreate the spin-polarized 
experiments of Loth et al. on single Mn and Fe atomtPSl 
Finally we test how the inclusion of the real part of the 
scattering self-energy modifies the SF-IETS for a non- 
polarized tip probing Mn monomers and trimers. 



II. THEORETICAL MODELS 

A. Interaction Hamiltonian and contour Green's 
functions 

We consider here the same sin gle- or bital tight-binding 
model used in our previous worka^SHZI to describe a mag- 
netic system (S) coupled to two non-interacting elec- 
trodes representing respectively the STM tip (tip) and 
the substrate (sub). The scattering region containing 
the magnetic nanostructure consists of a one-dimensional 
chain of N magnetic atoms. Each of the i-th atoms car- 
ries a quantum mechanical spin Si and it is character- 
ized by an on-site energy Eq. We assume that the tip 
and the substrate can only couple to one atom at a time 
in the scattering region thus to broaden the electronic 
level £o through the interaction with the electrode by 
rtip-s/sub-s- The Hamiltonian in the scattering region 
is then described by Hs — He + Hsp + i?e-sp, where Hg 
is the tight-binding electronic part, iJgp is the spin part 
and -ffc-sp describes the electron-spin interaction. The 
various terms can be written explicitly as 



Hc=eoY^ 4„CAa , (1) 

A a 

Hsp = Jdd 2_^ ^^ ' Sa+i+ (2) 

A 

N 

+ ^ {.g/iBB . S, + D{Si)' + EiiS'tr - (31)']} , 

A 

-f^e-sp — Jsd 2^ ic^^[cr]aa'C\a') ' Sa (3) 

A a, a' 

Aq 

where the electron ladder operator c}^ {cia) creates (an- 
nihilates) an electron at site i with spin a (a =tj i) ^nd 
on-site energy Eq- 

We model the spin-spin interaction between the local- 
ized {Si} spins by a nearest neighbour Heisenberg Hamil- 
tonian with coupling strength Jdd. Furthermore we in- 
clude interaction with an external magnetic field B (/^b 
is the Bohr magneton and g is the gyromagnetic ratio) 
and both uni-axial and t ransv erse anisotropy of magni- 
tude D and E respectively^MI, xhe electron-spin interac- 
tion Hamiltonian is constructed within the s-d modeP^^^Sl 
where the transport electron, s, are locally exchanged 
coupled to quantum spins, {Si} {d indicates that the lo- 
cal moments originating from the atomic d shell). In 
equation (|3| the interaction strength is Jsd and cr is a 
vector of Pauli matrices. The second term in Eq. (Is]) 
represents the potential scattering elastic contribution to 
the s-d interaction given by the exchange parameter Jq 
(note that this enters as a shift of the on-site potential of 
a given a tom) . The ratio, x = "^Z-^dj is typically in the 
range l-gHEU^ This term was not included in our pre- 



vious workalSITZl as it only becomes important for spin- 
polarized electrodes. 

In order to construct an electron-spin interacting self- 
energy we must first consider the KeldyslF^ contour- 
ordered single-body Green's functions (propagators) for 
both the electronic (G) and the spin (D) sub-systems in 
the electron-spin many-body ground state |) 

[G(r,r')]..'=-^(|TcK(r)4,(r')}|), (4) 

[Dir,r')U, = -z{\Tc{d^{r)dlir')}\). (5) 

These propagators describe a non-equilibrium system at 
zero-temperature. Here d}^ is a quasi-particle creation 
operator defined by the relation 



S\T)=ySl„^dl{T)d,,iT), 



(6) 



where i — {x, y, z} and m, n are eigenstates of iJgp of 
energy Smn- The matrix elements S'^„ = (m|S'*|n) de- 
termine the transition rates from the initial state n to 
the final state m. The quasi-particle operators are as- 
sumed fermionic in nature^''. Therefore the equilibrium 
population of a given state m is given by 

r^^ (7) 



Prr 



rl^ rl 
m ^m 



1 



where E is the energy, T the temperature and fee the 
Boltzmann constant. 

In the non-interacting case (Jgd = 0) the electronic 
system is not in equilibrium as the interaction with the 
electrodes establishes a steady state current. In contrast 
the spin system is in thermal equilibrium at the tempera- 
ture T. In this case the energy resolved lesser and greater 
Green's functions take the form 



[Gt{E)] 



[D>{E)Un = 



Pt=Ip-s(^)]- 



[S,^ 



sub— S 



{E)U' 



[nf(^)],„„ + [nf(s)]™ 



[E - Br. 



(kBTy 



(8) 



(9) 



In the electronic non-interacting Green's function [see 
Eq. ([8|] the coupHng of the sample to the tip and 
the substrate causes a broadening of the bare on-site 
level So of magnitude T = J2r,aa'[^rj-s]acy', where 
77 = {tip, sub}. The non-interacting self-energies take 
the form [^>iE)]„a' - [1 - fr,iE,V)][r^-sU' and 
[^<{E)],a' = UE,V)[rr,-sU', where ME,V) is the 
Fermi function in each of the 77-th leads at a bias V. 

In contrast the local spin Green's function of equa- 
tion (|9| describes a system, which is adiabatically cou- 
pled to a heat-bath of temperature T. This provides 
a very weak broadening of the single spin states Em 
of magnitude k^T. Such a heat bath keeps the spin- 
system in equilibrium and in the non-interacting case 
the population then resides mostly in the ground state 
[see Eq. Q]. For a ground state population of the spin 
system P^we have [U>{E)]„,n - <5,„„(1 - P°,)kBT and 



[n^iE)]r 



''mnPjn"'B-' ■ 



B. Spin-polarized electron self-energy 

In order to evaluate the effects that the interaction has 
on the electronic motion we must calculate the electron- 
spin self-energy. Here we take a perturbative approach 
and formally expand equation (HI) up to the n-th order 
in the interaction Hamiltonian, _ffp_=n, as 



[G{t,t')U 
(0|Tc{JTo-sp(Ti) 



y^HT^ 



n 



dri . 



dTr, 



c c 

.ffe-sp(T„)c.(T)cI,(r')}|0) 



, (10) 

[/(— 00, —00) 

where U is the time-evolution unitary operator and 
the time-averages are performed over the known non- 
interacting (Jsd = 0) ground state |0). The time inte- 
gration over r is ordered on the contour C going from 
—00 to -l-oo and then returning from -l-oo to — cx), since 
the ground state of the non-equilibrium system can only 
be defined at — cxPSl. 

In the following we consider the tip to have a spin 
polarization 77 (—1 < 77 < 1). This is defined as the 
spin asymmetry in the electronic coupling between the 
tip and the sample. An such the spin resolved electronic 
broadening is given by [Ftip-sln = ^'^rtip-s and 

[rtip-s]u = 2^rtip-s, where Ftip-s is the non-spin- 
polarized broadening. The substrate is assumed to re- 
main non-magnetic. As a result [G^lE)]^^ ^ [Gg'(£')]^^ 
so that we now must retain the spin indexes when con- 
structing the self-energy. The self-energy for the majority 
(^) and minority (J,) spins writes respectively as 



- J.'d E[^?(^ ± ann)]u^n(l " P^)\S'r, 



(11) 



and 



^UE)] 



rnn 



a 



G>(i?±f2„„)U;X 



I "nmX-' n'^mn ~r -Tn^i -'mjl'^mnl 



-J, 



2 Y^ 



sd / . 



[G>{E±n,r.n)]nPn{l~P.n)\S, 



mn I 



(12) 



The lesser (greater) self-energy describe an incoming 
(outgoing) electron that excites (relaxes) the spin system 
by ^mn , with a probability that depends on the occupa- 
tion of the spin levels P„j and P„ and on the spin selection 
rules 5*^+'- (note S+ = S"" + iSv and S- = S"" - iSy). 
The first term in both the equations (11 1 and ( [T2| ), pro- 
portional to Snrm Corresponds to the magnetoresistive 
elastic term of the s-d Hamiltonian of equation (l3|. The 
remaining contributions are inelastic in nature and de- 
pend on the spin orientation of the electron transferred 
from the tip. 



C. Spin-polarized spin self-energy 

When is magnetic and the tunnehng current tip car- 
ries a finite spin-polarization the spin system can be 
dragged out of equihbrium, in particular if the current 
density is intense. This mea ns th at the equilibrium con- 
ditions employed previousljff^'^, namely Pm=GS = 1 
and Pm^GS = is no longer valid. As a consequence 
we must now derive also an expression for the propa- 
gator and thus for the self-energy associated to the lo- 
cal spins. The derivation, up to second order in the 
electron-spin interaction is described in details in the 
appendix. In particular the total spin-self-energy also 
includes a zeroth-order contribution, which accounts for 
the non-interacting ( Jgd = 0) case. This is approximated 
by [U<{uj)U = P^knT and [Il>{u;)U = (1 - P^)kBT 
where P^ is the ground state population. Therefore, in 
absence of inelastic scattering, the spin system will re- 
main in thermal equilibrium with the heat bath and only 
the ground state will be occupied. 

By combining the zeroth and second order contribu- 
tions to self-energy we can write down a master equation, 
describing the non-equilibrium spin-population, in terms 
of the total self energy n^(^)^ 



dPn 
dt 



+ CX3 



IY. f dE{[Ii>{E)U^[D<{E)U„ 

m 

— oo 

-n<{E)Um[D^{E)]mn}- (13) 

After some rearrangement this can be written in more 
compact form as 



dPn 
dt 



Y. [Pnil - Pl)Wln - Plil - Pn)Wn 



{P° - P„)kBT , 



(14) 



where the bias dependent transition rate from an initial 
state |/) to a final state \n) is calculated after evaluating 
the integral described in equation (27) of the appendix. 
This finally writes 



Wnl = ~4 



(pJsd) 



{ 



vn' 



C,{Hn - Mr;' +r2/„)x 



'\Sni\ [r>)]n[rr,']n + l'5'n/l [r»)]u[r»;']u+ 

- I^n/nr„]u[r,']n + \sM^vM^v'\n} , (is) 



where C,{x) = x/ {l — e~^/'=BT~j ^^^ ^^ jg ^j^g chemical 
potential in lead 77 — {tip, sub}. Note that ({x) is such 
that for ?7 = 77' the resulting transition rates Wni are 
bias independent and do not contribute to the current. 
However, they do contribute to the spin relaxation time 
i.e. to the time taken by the localized spin system to relax 
back to its equilibrium state. Such a relaxation time is 



reduced if the coupling between the sample and the leads 
is increased. Furthermore, the smaller is the inelastic 
energy transition ^mn, the longer the spins will remain 
in the excited state before relaxing back to equilibrium. 
Finally, we note that the above expression is based on the 
assumption that the on-site energy is large enough for the 
density of states of the spin system to remain constant 
in the small energy window of interest. Therefore p = 
r/iel + T^). 



Returning to equation ( 14 1 note that we are only in 



terested in the steady state non-equilibrium population 
of the spin states at a given bias. Therefore we can set 
dPn[t)/dt = and reduce Eq. (14) to system of linear 



equations, which can be solved self-consistently. For an 
initial guess of the populations {Pi = P/*) we can it- 
erate Eq. (14), which define P„, with equation (27) of 
the appendix, which define the self-energy n> , until self- 
consistency is reached. We can then combine the re- 
sulting non-equilibrium population with the second order 
electronic self-energy calculated in Eq. (B to obtain the 
current. 



D. Real part of the electron self-energy 

In order to provide an explanation to the inherent 
asymmetry that has been observed in most of the STM 
experiments on magnetic atoms using a non-magnetic tip 
we return to the expression for the full retarded self- 
energy. This is defined by the Hilbert transform (note 
we will only consider the non-spin polarized case for sim- 
plicity) 



+00 

Wi^)^pv/f 



dE'Y.UE') + '>^UE') 



E~E' 



(16) 



-{S>,(i?') + S<t(i?')} 



By using the expressions derived in Section II C for the 
2"^^ order lesser and greater self-energies we can find an 
analytic expression for the real contribution to the re- 
tarded self energy 

Re[Si„t(i?)('^] - 2pJ^^ Y. I^™J'^"(1 - Prn)x (17) 

{kBTf- 



l{2.eo+Er.-sln[('' + ''™"^^" 



(P - fi™„ - /i^)2 -K {kBTf 



Such a final expression is heavily dependent on the on- 
site energy Eq but is also an odd function of the energy 
and the bias via its logarithmic dependence on the spin 
level eigenvalues with opposite polarity for -\-0.mn and 
— rim„. We will show in the results section that this is 
at the origin of the conductance asymmetry found in ex- 
periments. 



E. Electronic current 

We can finally unveil the effects that a non-equilibrium 
spin-population bares on the conductance profile of a 
magnetic nanostructure by taking the derivative of the 
current with respect to the bias voltage V. The current, 
Irj , flowing at the electrode ?/ 
as 



S=0 S=l 



S=2 



S=3 



{tip, sub} can be written 



Iv = 



+ CXD 



dEIr,{E) 



(18) 



T,{E) = |Tr{[I]> (i?)G<(i?)] - [T.<{E)G>{E)]} , (19) 

where G> (E) are the full many-body lesser/greater elec- 
tronic Green's functions. These are finally defined as 



[G^m 



[E^eo- Re[Si„t(£;)])2 -I- (r - Im[Ei„t(£;)])2 ' 

(20) 



and [J:tiE)] = [s|t(i?)](2) + [E|,(i?)](3) where the third 
order self-energies are calculated following Ref. [T7] . Here 
for simplicity we take the expression for the third order 
contribution to Sjnt obtained by neglecting any explicit 
spin-polarization. Such approximation is justified by the 
fact that the effects due to spin-polarization are small at 
the third order and that in doing so we avoid a rather 
cumbersome formulation. We do however consider the 
combination of non-spin polarized 3rd order effects with 
2nd order spin-polarized self energies to highlight subtle 
differences in the spectra. 



III. RESULTS 
A. Non-equilibrium population 

We start our analysis by first looking at the effects orig- 
inating from driving the spin system out of equilibrium 
with an electronic current. This attempts at explaining 
the experiments reported in Ref. 0], in which a STM 
tip (non-magnetic) is positioned above a Mn dimer de- 
posited onto a CuN substrate. The conductance spectra 
are measured for different tip to sample distances. Vary- 
ing the STM tip height is equivalent to changing both the 
current density and the electronic coupling between the 
tip and the sample. Non-equilibrium effects then appear 
as variations of the conductance profiles as a function of 
the STM tip height. 

Many of the parameters needed by our model can be 
extrapolated from a similar experiment carried on over 
Mn linear atomic chains deposited on CuN''. The five 
unpaired electrons in the Mn 3d shell suggest a 5 = 5/2 
ground state, as confirmed both by experiments and the- 
0]--^ |i4|24 | rp-j^g spin-spin exchange interaction between two 
Mn atoms is antiferromagnetic and has an estimated 
value of Jdd = 6.2 meV. As a result the ground state 
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FIG. 1: Excitation spectrum for an antiferromagnetically ex- 
changed coupled Mn dimer deposited on CuN. The ground 
state is a 5* = spin singlet. The first three excited state 
multiplets have respectively spin S — 1, S = 2 and S — 3. 
In the figure we also indicate the energy separation between 
the various spin multiplets in units of the exchange parameter 

Jdd- 



of the dimer is a singlet (total spin S — Q). The first 
excited state is a triplet with total spin S* = 1 and the 
energy splitting between the ground state and such first 
excited state is exactly Jdd- The next excited level is 
the quintuplet with total spin S — 2 and it is sepa- 
rated from the first excited state by 2 Jdd- This pattern 
continues throughout the spin manifold (see figure II]). 
The axial and transverse anisotropics are found to be 
D = —0.037 meV and E = 0.007 meV respectively. 
These cause the lifting of the spin multiplets degeneracy. 
The temperature is set at T = 0.5 K. The value of Jsd is 
estimated from density functional theory (DFT) to be of 
the order of 500 meWSl^ while Fgub-s is also found from 
DFT to be approximately 100 meV. In contrast Ftip_s 
remains an adjustable parameter with the chosen values 
ranging from 0.125 meV to 200 meV. Finally, in order 
to ensure a nearly constant density of states around the 
Fermi energy (Ep = 0) we set the on-site energy of the 
atom under the STM tip to be Eq = 1 eV. 

Figure [2] shows the conductance spectra obtained by 
simply taking the numerical derivative of the current 
[Eq. p^ ] with respect to the bias. We consider three dif- 
ferent tip to sample distance, corresponding respectively 
to weak (Ftip-s = 0.5 meV), intermediate (Ftip-s — 
50 meV) and strong (Ftip-s — 200 meV) electronic cou- 
pling. The evolution of the conductance lineshape as 
a function of Ftip_s is a direct consequence of the spin 
system being driven out of equilibrium. For Ftip_s = 
0.5 meV the STM tip is far enough from the sample 
to ensure that the spin system is always in its ground 
state between two subsequent electron tunneling events. 
Therefore the only transition detected in the G{V) pro- 
file is that between the S = ground state and the first 
excited state with S = 1. This has an excitation energy 
equal to Jdd and it does manifest itself as a conductance 
step at a voltage V — Jdd/e, with e being the electron 
charge. 

As the tip is brought closer to the sample (Ftip-s = 
50 meV) the first excited triplet level {S = 1) starts 
to populate. Now an incoming electron with sufficiently 
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FIG. 2: (Color online) Normalized conductance spectra for 
the Mn dimer calculated at different tip to sample distance, 
i.e. for different Ftip-s coupling strengths. We notice that the 
stronger is the coupling, the more the system is driven out of 
equilibrium. This results in the appearance of additional spin 
transitions, which manifest themselves as steps or drops of 
the conductance as a function of bias. 



large energy (2Jdd) can induce a second transition from 
from the first to the second excited state. Note that the 
S — 2 state is not accessible with a single electron tunnel- 
ing process from the ground state and it can be reached 
only if the spin system does not have enough time be- 
tween tunneling events to relax back to the ground state. 
For this case the transition appears as a reduction of the 
conductance at the critical voltage V = 2Jdd/e. The 
same spectroscopical feature is further enhanced at an 
even larger current density (Ftip-s = 200 meV), when a 
third conductance step appears at 3Jdd/e. This is asso- 
ciated to a transition from the 5' = 2 to the 5* = 3 spin 
state and it becomes possible only if the occupation of 
the 5 = 2 level is not zero, i.e. if the system is driven to 
this highly excited state. These results are in almost per- 
fect quantitative agreement with the experimental data 
(see Fig. 2of Ref. [4]). 

The evolution of the population of the various spin 
states (up to 5 = 2) as a function of bias is presented 
in figure [3J This is calculated in the case of strong tip 
to sample electronic coupling Ftip-s = 200 meV. In the 
figure one can note the strong spin-pumping from the 
ground state into both the first and the second excited 
state. The excitation to the 3'''^ excited state occurs at 
approximately 18 meV but is too weak to be observed on 
this scale. 



B. Spin-Polarized tip 

We now move on to consider the situation where the 
tip is magnetic, i.e. when the injected current is spin- 



FIG. 3: (Color online) Non-equilibrium population of the Mn 
dimer singlet (5 = 0), triplet [S — 1) and quintuplet {S = 2) 
states. The inset shows a magnified view of the population of 
the S — 2 state as it start to get populated at approximately 
12 mV. 



polarized. A gain we use as guide the experimental work 
of Loth et aZ.SEl. The STM tip is now spin-polarized by 
placing an additional Mn atom at its apex, while also ap- 
plying a strong magnetic field perpendicular to the sub- 
strate of 3 T. In this case the spectrum is collected from 
a single Mn or Fe ion on the surface (not from a dimer). 
In the case of Mn, the atom exhibits a weak anisotropy 
on CuN (see previous section) and the strong magnetic 
field effectively produces a Zeeman split of the six levels 
of the S — 5/2 Mn spin manifold. The direction of the 
magnetic field in these experiments is chosen so that the 
ground state of the Mn spin corresponds to the magnetic 
quantum number m — -H5/2. Since the same magnetic 
field is applied to the Mn atom on the tip's apex, the tip 
and atom are both spin-polarized and coUinear. 

Figure |4] shows the calculated spectra for the system 
described above. In particular we consider magnetic field 
strength of 3 T and either weak (Ftip_s = 0-5 meV), 
intermediate (Ftip_s = 50 meV) or strong (Ftip-s — 
200 meV) tip to sample couplings. The on-site energy 
is fixed at e^ — 2eV and the value of Jgd = 500meV is 
infered from the work of Lucignano et ai^^. The tip 
spin-polarization constant and the inelastic ratio that 
best fit the experimental data are respectively 77 = —0.3 
and X = 1-5- In the weak coupling regime (when the 
local spin remains always close to equilibrium) the lo- 
cal spin resides almost entirely in its ra = +5/2 ground 
state. Due to the spin-exchange selection rules and to 
the coUinearity of the tip and the sample, only the mi- 
nority carriers can excite the local spin out of the ground 
state. For a tip spin-polarization of 77 = —0.3, there are 
more minority electrons coming from the tip than those 
coming from the substrate. As a result, the intensity 
of the inelastic interaction will change depending on the 
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FIG. 4: (Color online) Normalized conductance spectra for a 
single Mn atom explored by STM at different tip to sample 
electronic couplings Ftip-s- Here a magnetic field of 3 T is 
applied along the z-axis. The asymmetry in the conductance 
profile is is due to the spin-polarization of the tip. Such an 
asymmetry is more pronounced as the system is driven further 
away from equilibrium. Comparison is made between second 
and third order calculations. 



direction of the current. This creates an asymmetry in 
the conductance spectrum with respect to the applied 
bias. The additional Hneshape features appearing in the 
weak coupling case (the conductance decay following a 
conductance step) are due to the third order Kondo-like 
self-energy, which produces a logarithmic decay at the 
conductance steps. This result is in good agreement with 
experiments (see Fig. 4 of Ref. H]). 

When the spin of the Mn ion is driven further out of 
equilibrium, in particular in the strong coupling case, the 
bias asymmetry becomes more pronounced, such spin- 
pumping phenomenon can be appreciated by looking at 
figure IIIB| [A), where we show the populations of the six 
spin states of the Mn atom as a function of bias for strong 
tip to sample coupling (Ftip-_s = 200 meV). From the 
figure one can see that as the bias increases the to — +h/2 
ground state gets depleted in favour of populating the 
other five excited states. In particular already at T^ ~ 
10 mV the population of the m = —5/2 level is larger 
than that of the ground state. In th figure we also plot 
the average magnetization, which is defined as {S^) — 
Ylm PmS^n^ [scc panel (b)]. Intriguingly we find that for 
negative biases the spin is effectively flipped from m — 
-1-5/2 to TO = —3/4 over 25meV range. Such spin flipping 
results in a large dip in the conductance for negative 
biases as the tip is no longer coUinear to the sample. 

In Fig |6] we present the calculated spectra for the Fe 
atom in the spin polarized case. We choose parameters 
in this case that conform with the experimental data of 
Loth et al. in Ref.!^. The Fe atom is assumed to carry 
a quantum mechanical spin of 5' = 2 and it also exhibits 
a transverse easy axis anisotopy of -D = — 1.53meV and 
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FIG. 5: (Color online) (a) Non-equilibrium population of the 
various spin states of Mn on CuN as a function of bias and (b) 
the resulting average magnetization. These have been calcu- 
late for a magnetic field of 7 T aligned along the z-axis and for 
strong tip to sample electronic coupling Ftip-s = 200 meV. 
We notice that as the spin is driven far away from its equilib- 
rium ground state the magnetization flips its direction. 



an axial anisotropy oi E = O.SlmeV. We again assume a 
large value of the onsite energy, 2eV, and we examine the 
spectra in the strong coupling case of Ftip_s = 200meV 
with a tip polarisation of ?7 = 0.35 and magnetic field 
strength 3T as used in experiments. In Fig 6(a) we 
present the conductance spectra for the two cases when 
the magnetic field is parallel and perpendicular to the 
easy axis of the atom (the z— axis in this model). As pre- 
viously, we present this for both second and third order 
calculations. Firstly, we notice that the spin polarized 
tip affects the spectra only in the case of parallel mag- 
netic field where a clear bias assymetry is produced. No 
significant assymetry is found in the perpendicular case. 
This conforms with the experimental findings and is due 
to the fact that electron spins in the tip are no longer 
colinear with localised spin of the Fe atom. 

As found in previous workal^ the inclusion of third or- 
der effects is vital in reproducing the corresct logarithmic 
decay at each of the conductance steps, which is par- 
ticularly noticable for the perpendicular magnetic field. 
More significantly, experimental spectra for the parallel 
case exhibit a zero bias conductance dip which is absent 
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FIG. 6: (Color online) (a) 2nd and 3rd order conductance spectra 
for the Fe atom with spin polarization r] = 0.35. The magnetic field 
of 3T is applied both parallel and perpendicular to the easy axis 
of the Fe atom, (b) Second dervative of the current for the same 
spectra shown in (a). 



We test this approach by considering the case of a non- 
spin-polarized tip and a single Mn atom. We use the same 
anisotropy parameters as for the Mn dimer but, for the 
sale of simplicity, we keep the spin always its equilibrium 
state and choose Ftip-s = 0.5meV. Figure M shows the 
resulting conductance spectra for three different choices 
of the on-site energy e. It is cear that the closer e is to the 
Fermi energy (0 eV), the greater is the bias asymmetry, 
while as e is increased, the conductance profile becomes 
more symmetric. In this respect, the formalism outlined 
here is in agreement with the Fano lineshape argumenP^l 
where the degree of asymmetry for electrons tunneling 
through a single impurity is given by a ratio of the real 
to the i mag inary contributions to the interacting Green's 
functioiPl 
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in the 2nd order spectra but appears strongly when third 
order effects are included. This can also be seen the cal- 
culation of the second derivative of the current in Fig 
6(b) where a clear zero bias anomaly is evident in the 
third order case. 



FIG. 7: (Color online) Normalised conductance spectra for a 
single Mn single calculated by including the real contribution 
to the interacting electron-spin self-energy for different on-site 
energies £o- Note that the conductance asymmetry increases 
with decreasing eo, i-e. as the onsite energy moves closer to 
the Fermi level. 



C. Non-spin polarized asymmetry 

We finally move to discuss the inherent asymmetry 
measured in the conductance profile, which is usually ob- 
served even if the tip is not spin-polarizecP^. We model 
this lineshape feature by including the real part of the full 
interacting electron-spin self-energy in the description 
[see equation (17)]. The structure of this contribution to 



the self-energy shows an explicit dependence on the on- 
site energy, Eq, and also a logarithmic peak of width fceT 
at the onset of an inelastic transition {E — ^^i = i^mn)- 
The asymmetry arises from the difference in polarity of 
the logarithmic peak for ±ri„i„. The self-energy is thus 
an odd function of both energy and bias. This results 
in the conductance profile having a decrease of the step 
heights for V = —Q„in/s and a increase of them for 

V = +nmn/e. 



As an additional test we consider the case of a Mn 
trimer, whose spectrum was shown first by Hijibehedin 
et al™ to exhibit a large bias asymmetry when mea- 
sured with a non-magnetic tip. We model this system 
by choosing an antiferromagnctic nearest neighbour ex- 
change coupling J^J =2.3 meV. Furthermore, in order 
to accurately describe the position of the principle con- 
ductance steps in the conductance profile, we also include 
a ferromagnetic second-nearest-neighbour interaction be- 
tween the local spins of magnitude J^J =-1.0 meV^. 

We again choose to keep the spin system in equi- 
librium and therefore consider weak coupling between 
the STM tip and the second atom in the trimer chain 
(ftip-s — 0.5 meV). The best fit to the experimental data 
is found with Eq — — leV. Figure ?? shows the model 
fit to the experimental data (from Ref. [3j)). Whereas 
previous calculations did not predict any conductance 
asymmetrylSl it is clear from the figure that the inclusion 
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of the real part of the self-energy in the description pro- 
duces a significant conductance asymmetry. This is most 
prominent at the principle step height (« ±16.5nieV) for 
each bias polarity. Although the step height for the nega- 
tive bias is not as small as that found experimentally, the 
qualitative trends are similar. In particular we notice the 
logarithmic conductance increase (reduction) that occurs 
before (after) the onset of the step at F = +16.5meV, 
which also originates from the third order contribution 
to the self-energy. 

In this work, based on a perturbative approach of the 
s-d model, we have shown that the entire lineshape de- 
scription can be re-conciliated with experiments by con- 
sidering an expansion of the self-energy to the third order, 
which also includes its real part. As such we have shown 
that the conductance asymmetry can be described also if 
the electronic orbitals forming the samples spin are not 
explicitly taken into account as also suggested by Del- 
gado and Fernandez-Rossiep22l. 




-10 10 20 
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FIG. 8: (Color online) Comparison between the experimen- 
tal (red) and theoretical (black) conductance spectra for a 
Mn trimer on CuN probed by a non-magnetic STM tip. In 
the calculation we have included the real part of the interact- 
ing electron-spin self-energy. This provides the conductance 
asymmetry with bias. 



IV. CONCLUSIONS 

In conclusion we have studied the lineshape details of 
the conductance profile of Mn atoms deposited on CuN 
and probed with a STM tip either or not carrying spin- 
polarization. In particular we have looked closely at the 
asymmetry of the conductance with the bias polarity. 
Firstly, we have extended our perturbative approach to 
spin-scattering to the spin-polarized case and considered 
an expansion of the complex part of the electronic prop- 
agator up to the third order. This allows us to reproduce 
the logarithmic decay of the conductance subsequent a 



conductance step, which is observed in experiments but 
could not be explained by a second order theory. 

When the current density is increased and the tip is 
spin-polarized the conductance profile starts to develop 
a significant asymmetry with respect to the bias polar- 
ity. These are indicative of the spin system being driven 
out of equilibrium. We have then derived a second order 
expansion of the spin-propagator capable of evaluating 
the non-equilibrium population of the various spin energy 
levels. This was put favorably to the test against a series 
of experiments probing a single Mn and Fe ions with a 
spin-polarized STM tip in an intense magnetic field. Fur- 
thermore the same formalism was capable of describing 
excitations occurring away from the ground state for a 
Mn dimer probed by a non-magnetic tip. Also in this 
case the agreement with experiments is very satisfactory. 

Finally, in an attempt to describe the bias asymme- 
try in the case of non-spin-polarized STM tips we have 
derived an analytic expression for the real part of the 
electron-spin interacting self-energy. This contains loga- 
rithmic peaks at the excitation energies that are odd with 
respect to energy and voltage. Such parity results in an 
asymmetry in the conductance profiles. Such a scheme 
was tested for the case of a Mn monomer and a Mn trimer 
and compares reasonably well with experiments. 
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VI. APPENDIX 

Here we wish to derive a method for calculating the 
steady-state non-equilibrium distribution of the spin en- 
ergy levels populations, Pm{V), due to the coupfing with 
the electrodes. In order to do so we expand equation (Is]) 
up to the n-th order in the interaction Hamiltonian 



[D{T,r')U,^Y. 



(-0 



n+l 



dri 



dTn X 



c c 

(0|rc{ge-sp(ri) . . . ge-sp(T-„)rf„(T)dt^(TO}|0) 

C/(— oo, — oo) 



, (21) 



where U is the time-evolution unitary operator and the 
time-averages are over the known non- interacting {Jsd — 
0) ground state |0) . As in equation ( 10 ) the time integra- 
tion over r is ordered on the contour C going from — oo 
to -|-oo and then returning from -|-oo to — ocPS. 
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By inserting the expression for Hc-sp from equation 
([3]) into the equation above and by expanding up to the 
second order we obtain [note for the ease of the descrip- 
tion we omit the elastic contribution of Jq, which is then 



included in the final expression in equation (26)] 



[D(7 



r'^l(2) = 



2! 



-J: 



sd 



E 



Q,a' ,/3,/3' ^ 



dri / dT2 X 



irc{4(ri)c„,(ri)4(r2)c^,(r2)d„(r)dt„(r')}|0) 
J2{0\Tc{S\n)S^iT,)}\0)[a\^,[a'W. 



(22) 



where the indices i and j run over the cartesian coordi- 
nates X, y and z for the given spin coupled to the tip (the 
tip make electronic contact with one spin only). We now 
substitute into equation (p2| ) the operator breakdown of 
the spin from equation (161 



[D{, 



/Nl(2) _ (_-*)^ 72 



2! -^^d 2^ 



dri / dT2 
c 



x(0|Tc{d„(T)4(Tl)4'(Tl)4(T2)dp(T2)dL(^')}|0) 

X Y. (0|Tc{c„,(Ti)4(r2)c^.(T2)ct,(ri)}|0) 



J2si,,siAc^^U'W'] 



/3/3' 



(23) 



The time-ordered contractions of the two brackets in 



equation (23) can be re- written in terms of their re- 



spective non-interacting Green's functions, Dq(t,t') and 
Go(t, r') as follows 



-•^s^d E J ' 



\D(T^T')\y^l,^^Jt^ ^ /dn \dT2 

' '' c c 

X <5nfe^i/c"^mr[-Do(T, ri)]„„[i:»o(Ti,T2)]/;[I?o(T'2,T')]mm 
X 2_^ 5a'P^aP'\GQ{Tx,T2y\fifi\G{){T2-,Ti)\aa 



J2si,,siA'y'U'[T' 



1/3/3'' 



(24) 



where the extra factor of 2 emerges from the fact that 
a second contraction of the time-ordered bracket merely 
exchanges ti and T2. Then, by using Dyson's equatioiP^, 
one can write the second order contribution to the inter- 



acting spin self-energy (11). This reads 

[n(ri,T2)]i2)^ = -2j2^^[Go(Ti,r2)]^;5[Go(r2,n)], 



a,l3 



X 5][Do(Tl,r2)]z,^^5;^5LK]„^K 



|/3a 



(25) 



where we have evoked the assumption that the elec- 
trons are spin degenerate thus omitting the spin index 
on Go('''i,'''2) and including a factor of 2. We now calcu- 
late the real-time quantities, such as the lesser (greater) 
self-energies, by using Langreth's theorem for the time 
ordering over the defined contouifSSl. After including the 
elastic contribution we obtain 



[n^(ti,t2) 



(2) 

nm 



a,P 

I 

{siiSlJa^U[a']po. + 5,,5^pxSlmW%p) ■ (26) 



i,j 



On computing the Fourier transform we note the two 
different expressions for the lesser and greater Green's 
functions are 

+CXJ 

^ j du[G<{u:)]p^[G>{Lo±{E~ei))Uy. 

"./3-00 

E ('5'ni'S'LK]a0[o-^]/3a + 5»j <5a/3X'5'^m K]a/3) , (27) 

where we have defined P^ — Pi and Pf" = 1 — P; and 

P>o{ti,t2)ii = P;>exp[-je;(ii -^2)/^]- By assuming that 
the spin system is in thermal contact with a heat bath 
kept at the temperature T, the energy levels ei should 
be broadened by the factor /? = fceT. This can be ne- 
glected for the ease of the calculation since in general 
T << 1. However we do not disregard the broadening 
in the electronic Green's function due to contact to tip 
and substrate as this is pivotal to the calculation of the 
non-equilibrium spin populations. 
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